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Abstract. We discuss mechanisms for angular momen- 
tum transport in the clumpy medium of the circumnu- 
clcar disk at the Galactic center. The viscosity due to 
clump-clump collisions is found to be less than some crit- 
ical viscosity; this meets the conditions at which a collec- 
tive mode of nonaxisymmetric shear perturbations in the 
disk is able to grow until going into the saturation regime 
where fully developed turbulence is established. We hnd 
that the angular momentum transfer due to this turbulent 
viscosity turns out to comparable to the transport due to 
magnetic torques. Taken together, the turbulent and mag- 
netic transfer of angular momentum are able to provide 
the inflow of mass into the central parsec with a rate of 
about 10 -2 M yr~ 1 , consistent with the available data. 

Key words: Galaxy: center of - 
- interstellar medium: clouds 



accretion, accretion disks 



10 -2 M /yr on the scale of hundreds pc. Earlier, a similar 
mass inflow rate was inferred to exist on the scale of the 
circumnuclear disk, i.e. at 1.5 pc <i?<10 pc (Genzel & 
Townes 1987, Jackson et al. 1993). 

We also note that the star formation in the inner region 
of the Galaxy requires such a mass supply. 

However, appropriate physical mechanism(s) which 
would be able to provide the necessary rate of momentum 
transfer within the central 10 pc or so is (are) obscure 
so far. The present paper aims at consideration of such 
mechanisms. In Sec. 2 and Appendix, we discuss viscos- 
ity in a clumpy disk due to clump interactions, both with 
and without self-gravity included. Sec. 3 deals with short- 
wave instability of the clumpy, viscous disk; the necessary 
condition for developing small-scale turbulence is estab- 
lished here. In Sec. 4, we apply this mechanism to the 
circum-nuclear ring at the Galactic center. The results of 
the paper are discussed in Sec. 5. 



1. Introduction 

Transport of angular momentum in the disks of spiral 
galaxies is one of the central issues in galactic dynam- 
ics. It is also of prime importance for fueling the central 
engines in active galactic nuclei (AGN). With regard to 
our own Galaxy this topic has been recently explored in a 
scries of papers (von Linden et al. 1993a, b; Biermann et 
al. 1993) with the conclusion that the energy input from 
the SN explosions can feed the turbulence of the interstel- 
lar matter so as to provide the effective viscosity high as 
required to feed the star formation. 

An immediate implication is that the viscous trans- 
port of angular momentum outwards is accompanied by 
inflow of gas inward with the rate that might be as high as 

Send offprint requests to: L. Ozernoy 



2. Viscosity due to cloud-cloud interactions 

In a recent version of the unified model for fueling AGN, 
Begelman et al. (1989) proposed that the inflow of mat- 
ter (driven by global axisymmetric gravitational instabil- 
ities on large scales) proceeds on small scales in the form 
of a 'disk' composed of randomly moving clouds, which 
are embedded in a low-density medium with a small fill- 
ing factor. There is an essential uncertainty in our knowl- 
edge about the overall configuration, dynamics, and con- 
finement mechanisms for thermal gas clouds observed in 
AGN, but as for our Galaxy, such a disk is known to exist 
as the circum-nuclear disk, or ring that is rather clumpy, 
indeed (Jackson et al. 1993). Begelman et al. (1989) as- 
sumed that the viscosity in the disk is provided by colli- 
sions between clouds. Below, this collisional mechanism of 
viscosity is compared with some others. 
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We consider the following simple model of a cloudy 
disk: The clouds are orbiting in an external gravitational 
field and have some random peculiar velocities. It is as- 
sumed that the clouds have a small filling factor and are 
embedded in a low-density medium that provides a con- 
finement of the clouds. In the differentially-rotating cloudy 
disk, the angular momentum is transported due to cloud- 
cloud interactions (which include, and are not just re- 
stricted to, collisions). Let us address the shear viscosity 
associated with these interactions. 

The most elaborated models to calculate the shear 
viscosity have been considered by Goldreich & Tremaine 
(1978) and Stewart & Kaula (1980), hereinafter referred 
to as GT and SK, correspondingly. GT considered con- 
tact inelastic collisions of non-gravitating spherical parti- 
cles obeying an anisotropic distribution function, whereas 
SK considered gravitational (elastic) encounters of parti- 
cles obeying a Maxwellian distribution. In the both cases, 
as shown in Appendix A, the viscosity coefficient can be 
represented in the form 



mean number of interactions suffered by a cloud in passing 
through the disk. More accurately, 



t ~ Tra 2 nh 



(2.5) 



assuming a to be the largest of geometrical and gravita- 
tional influence sizes. Here h is the thickness of the disk 
given by 



a v 



(2.6) 



Since h ~ S/mn, where E is the surface density of the 
disk, Eq. (2.2) can be rewritten as 



r ~ 



na 



(2.7) 



which implies one more interpretation for r: it is the cov- 
ering factor, C, or the fraction of disk area covered by 
clouds when they are placed as a monolayer. 

The filling factor of the system of clouds, i. e. the frac- 
tional volume filled by the clouds is 



A, 



n b 



2_2 



1 



(2.1) 



Here a v is the one-dimensional velocity dispersion, fi is the 
orbital angular velocity in the disk, r = Q.U is the 'optical 
depth' to cloud-cloud interactions (see below), U is the 
free path time (index i = c, G; c stands for collisional, 
or contact, interactions, G stands for gravitational ones), 
and Ai, Bi are the constant coefficients defined below. 

Collisions between the clouds result in diminishing a v , 
whereas gravitational encounters tend to increase it. For 
cloud-cloud collisions, 



t c = ("7ra 2 ner„) 



(2.2) 



where a and n are the typical size of a cloud and the spa- 
tial number density of the clouds, respectively. For gravi- 
tational encounters between the clouds 



to = 



4^/irG 2 m 2 n 



(2.3) 



(Braginskii 1965), where m is the typical mass of a cloud. 
Evidently, ac — 2 (3t/tt) Gm/a 2 can be considered as 
an effective size of the domain for gravitational influence 
of the cloud. 

Coefficients Ai and Bi in Eq. (2.1) take the following 
values: 



A c 
A G 



0.46; 

{1.25 if v = const, 
1 



0.83 if Q. oc r- 3 / 2 



B c = 0.97 
Bg = 1.95. 



(2.4) 



The optical depth r of the disk is a convenient param- 
eter describing how effective are the interactions between 
the clouds. By order of magnitude, it is nothing but the 



F = -Tra 3 n . 



(2.8) 



Eqs. (2.5), (2.6), and (2.8) yield one more expression for 
t containing F: 



TTa 2 ncr v 3^/i 

t ~ = —b — 

fl 4 a 



(2.9) 



Evidently, F <C 1 for any cloudy disk with a <C h unless 
r'Zh/a > 1. 

Here we emphasize that cloud-cloud collisisons cannot 
be discussed without noting that the magnetic fields per- 
meate clouds, and are likely to make such collisions more 
efficient by increasing the effective cross section: When 
two clouds collide, it is unlikely that they just slide along 
a given straight flux tube. First of all, if this imagined 
flux tube were not exactly on a circle, then by angu- 
lar momentum conservation the clouds would not go in 
a straight line, and second, by virtue of the general dis- 
tribution of velocities it is very unlikely that two clouds 
would just match in proper velocities to be able to slide 
along a flux tube, and, third, the energy density in flux 
tubes is unlikely to sufficiently overpower the kinetic en- 
ergy of clouds to do this. It follows that it is indeed likely 
that the flux tube will be twisted, thus strengthened in 
magnetic field, and therefore the clouds may interact even 
at some distance. This means that cloud-cloud collisions 
may involve a larger effective cross section than just the 
geometry would imply. 

Furthermore, with magnetic loops and reconncction in 
the region above the disk, the effective scale height may 
well be larger than the scale height of the visible cloud 
distribution. 

It follows that the estimate above may be an underesti- 
mate just as well as an overestimate; the observational fit 
and interpretation given to the data by von Linden et al. 
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suggests that the viscosity derived above for cloud-cloud 
collisions is an underestimate. 

It is instructive to compare the viscosity coefficient in 
a cloudy disk [Eq. (2.1)] with that in a typical thin, but 
continuous disk: v ~ v s l, where v s and I are the sound 
velocity and the mean free path length, correspondingly 
Qualitatively, in a continuous disk I is anticipated to be 
much smaller than I in a cloudy disk. If v s ~ a v , the 
viscosity in a cloudy disk exceeds that in a continuous 
disk, bearing in mind some reasonable assumptions about 
the disk parameters. Before specifying them, we would like 
to discuss one more mechanism for viscosity in a cloudy 
disk, this time of a collective origin, proposed recently by 
Fridman & Ozernoy (1992). 

3. Small-scale turbulent viscosity 

We consider, for simplicity, both the typical size of the 
clouds, a, and the mean free path of the clouds, I, to obey 
an inequality: 

a, I < h . (3.1) 

The hydrodynamics of a cloudy accretion disk with respect 
to 2-D nonaxisymmetric shear perturbations is similar to 
that for a disk with the continuous distribution of the 
matter explored by Fridman (1989) who found the solution 
of this problem in a small-amplitude, short-wave limit 



(3.2) 



where k is the wave vector. In a local rotating coordinate 
system, the solution for the radial component of the per- 
turbed velocity, v\ r , reads: 



1 + /3 2 



Vlrjt) _ 

Mo) ~~ i + {p + tf 



exp(-i//z/ cr ) 



where 



n MO) - A 

(3=-^, t = At, A 



dr 



A[k 2 p t(l + /3 2 +l3i+-t 2 ) 



(3-3) 



(3.4) 



(3.5) 



Here f2(r) is the angular velocity of the disk, t is the time 
elapsed since the perturbations were "turned on" , and v CT 
is a (timc-dependent) critical viscosity. When v <C v CI , 
viscosity does not play any role in the disk dynamics. 

The shortwave perturbations under consideration be- 
have as incompressible modes. The solution (3.3) describes 
how the vorticity decays with time in a viscous fluid. Some 
limiting cases of interest can be revealed from Eqs. (3.3)- 
(3.5). As t — > oo, there appears the asymptotic solution 
~ exp(— at 3 ), which has been known earlier (e. g. Timo- 
feev 1976, Zaslavskii et al. 1982). Another limiting case 
when viscosity is absent {y = 0), deserves a more detailed 
consideration. 



3.1. Non-viscous case (v = 0) 

In the absence of viscosity, the solution (3.3) goes into 

v lr (0)(l + [3 2 ) 



Vlr{t) = 



! + (/? + t) s 



(3.6) 



which was obtained by Lominadze et al. (1988). In this 
case, solution (3.6) implies that 



vi r (t)k 2 j_(t) = v lr (0)k 2 j_(0) = const , 



(3.7) 



where 



fci(t) = k 2 r (t) + kl, k r (t) = k r (0) 



kpt . 



Eq. (3.7) is a 2-D analog of the Thomson theorem on the 
conservation of vorticity in an incompressible fluid flow 
(Fridman 1989). 

The solution (3.6) is shown in Fig. 1 by the solid lines. 
Evidently, the velocity perturbation is growing (implying 
instability) if the denominator in (3.6) is decreasing. At 
the moment when the denominator has a minimum, i.e. 
/? + t* = 0, where t* is given by 



k r (0) 



Vir(t) reaches its maximum equal to 



^lr max(^*) 



fc r (0) 

km 



Mo) 



(3.8) 



(3.9) 



and then goes to zero at t — > oo. The growth of v\ r at < 
t < —(3 is a result of a decrease of k±(t) while vi r (t)kj_(t) 
keeps constant due to the Thomson theorem (3.7). 

The minimum of the denominator in Eq. (3.6) at 
(J3 + U) = implies that {k r (0)/k v + r \dQ/dr\ t) = 
as d£L/dr < 0. Since r > 0, t > 0, the necessary condition 
for the perturbations to grow with time is given by 



fc r (0) 

kin 



< . 



(3.10) 



Therefore, in a differentially rotating disk, the growing 
short-scale spiral perturbations can only be leading. 

To clarify possible misunderstandings, we note several 
points in the following: 

In Appendix B, we show that in a differentially rotat- 
ing disk there are no incompressible short-wave perturba- 
tions oscillating with the epicyclic frequency at all; in a 
solid-body rotating disk such perturbations are stationary, 
in agreement with our eqs.(3.6)-(3.10). 

In Appendix C, we demonstrate that Coriolis forces 
are irrelevant to the solution for velocity perturbations in 
a plane shear layer, which we found to be similar to short- 
wave, low-frequency perturbations in a rotating gaseous 
disk. We emphasize that Coriolis forces have nothing to 
do with the physical meaning of our solution. The latter 
is associated with vorticity conservation (in the limit of 
small viscosity). 



4 



L. Ozernoy et al.: Angular momentum transport 



We consider perturbations with wavelengths smaller 
than the disk's width. Consequently our case is close to a 
cylinder fluid layer. Meanwhile the short-scale perturba- 
tions under the condition (3.2) are incompressible (Frid- 
man 1989). Therefore analytical results obtained in this 
paper as well as in our previous papers (Fridman 1989, 
Fridman and Ozernoy 1992) lend to support interpreta- 
tion of well-known laboratory experiments for a liquid 
flow between two rotating cylinders. The latter demon- 
strate a nonlinear instability and developed turbulence 
(e.g., Lukashchuk & Predtechenskii 1984). Therefore, we 
have all the reasons to believe that, for short-wave per- 
turbations, a nonlinear instability is likely to occur. Like- 
wise, in Appendix 3, we show that the dynamics of short- 
wave incompressible perturbations in a rotating disk is 
completely equivalent to that in a plane shear layer. This 
proves the equivalence of the nonlinear instability in a 
plane shear layer and in the situation investigated here. 

It is instructive to compare the above instability with 
that considered by Goldreich & Lynden-Bell (1965) (re- 
ferred to as GLB hereinafter) for a self-gravitating, dif- 
ferentially rotating disk. There are two basic differences 
between our and their situations. First, different branches 
of growing perturbations are considered: we deal with the 
vortex branch while GLB do with the sound-gravitational 
branch. Second, the conditions for the growth of perturba- 
tions are different, too. Under condition (3.10), the lead- 
ing vortex perturbations grow, due to the conservation of 
vorticity, in an otherwise stable disk. GLB deal with the 
disk stable for axially-symmetrical perturbations, and the 
growth of their non-axially-symmctrical, long-wave per- 
turbations proceeds due to a leading role of self-gravity. 
Toomre (1982) christened this growth as "swing amplifica- 
tion" . It is worth noting that, in our case, a similar swing 
amplification occurs for the (vortex) short-wave perturba- 
tions and it does not require self-gravity. 

3.2. General case (non-zero viscosity) 

The growth of perturbations with accounting for viscos- 
ity, which is described by Eq. (3.3), is shown in Fig. 1 by 
the dotted and dashed lines. The "instability" described 
by the exact solution (3.3) is of very specific kind as v\ r 
tends to zero at i — > oo. Nevertheless, during a finite 
time U given by Eq. (3.8) the perturbations arc in the 
growing regime and their amplitude increases by a fac- 
tor of k 2 ± {0)/kl(-P) = (J3 2 + 1), which can be > 1 if 

\/3\ = |Mo)/M » !• 

Let us suppose that by the moment t = t* = \f3\ 
when the amplitude of any perturbation arrives at its 
maximum, the viscosity does not play any essential role, 
i.e. v <C v CT \t^). Though this growth in the amplitude 
proceeds during the time interval At ~ only, it could 
onset the local turbulence. In this case, turbulence can be 
established everywhere in the disk as a superposition of 
the spiral perturbations originated in different points of 



the disk and on different moments of time; all of them 
have experienced a growth during At ~ t*. 

The turbulence appeared as a result of the growth of 
the perturbations can be characterized by some turbu- 
lent viscosity, which generally is much larger than the 
molecular one. It is possible that a steady-state regime 
will emerge in which 



V —* f C r(t*) = fturb 



(3.11) 



In this regime, the decay due to viscosity is strong enough 
to provide a steady-state level of the turbulence so that 
the amplitude of the perturbations is kept more or less 
constant in time, i. e. 



Vlrjt) 

Mo) 



1 at 



P(t.) = -/3~ 1 



(3.12) 



Substituting Eqs. (3.11) and (3.12) into Eq. (3.3) with 
taking into account Eq. (3.5) and assuming that, by order 
of magnitude, A ~ ft, we arrive at the following transcen- 
dental equation: 



exp -- 



4 fturb k 2 



ft 



Taking the logarithm of this equation one finds: 

ft 



^turb — 0.5 



U2 



(3.13) 



(3.14) 



It is straightforward to see that the basic contribution 
into the turbulent viscosity is given by the perturbations 
of the smallest fc^'s; therefore f tur b — 0.5ft/(fc^,) m j n . From 
Eqs. (3.2) and (3.12) one has \k v \ ~ \k r \ » ft, -1 , whence 



\k 



h- 



i. e. 



t'tm-b — 0.5 flh 2 



(3.15) 



By substituting h = a v /il into Eq. (3.15) the latter can 
be rewritten in the form: 



fturb ^ 0.5 



ft 



(3.16) 



It is instructive to compare the expression for the 
Bohm diffusion coefficient (Db) with our Eqs. (3.14) and 
(3.15) for ^turb- A well-known estimation of the Bohm dif- 
fusion coefficient for a strong turbulence plasma is given 
by (see e.g., Kadomtsev, 1964, or a review by Horton, 
1984): 



{D B \ 



(7l)e 



j_ / mm 



where is the maximum linear growth rate of the drift 
instability, and (fcj_) is the minimum wave number, also 
from a linear theory. They are much the same: in our 
case (7L)max ~ 0.5 ft and instead of k± we substitute 



~ h- 



Kadomtsev obtained (I?s) ma x by using 



the relationship k^p <~ 1, where p is the Larmor radius. At 
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k± ~ p^ 1 , Db has a maximum. But the Larmor radius in 
plasma corresponds to the epicyclic radius in graviphysics. 
The latter, in fact, is the thickness of disk, k± ~ h^ 1 , and 
this, actually, is being used here. 

This viscosity whose value is given by Eq. (3.15) or 
(3.16) was called by Fridman & Ozernoy (1992) anomalous 
in the same sense as one introduces anomalous resistivity 
and anomalous diffusion in plasma physics: The origin of 
this viscosity is in fully developed turbulence which is es- 
tablished in the saturation regime described above. 

Eq. (3.15) or (3.16) could be derived from dimensional 
arguments (of course, without the numerical coefficient) as 
an estimation of viscosity in a rotating disk with turbulent 
motions. However, we should emphasize that without an 
analysis such as one given above it would be impossible to 
reveal an underlying physical mechanism for the origin of 
such turbulence. 

It is important to test whether the necessary condition 
for developing of small-scale turbulence 



v < y c 



(3.17) 



is met. To this end, let us find the ratio v/vtuxb as a func- 
tion of t = Clti, i. e. of the basic parameter that charac- 
terizes the number of interactions per one revolution: 



v 

v„ 



2A t r 
BfT* + 1 



(3.18) 



Asymptotically, 



v 

Vex 



2A t T < 1 if T < 1 , 
^T" 1 < 1 if T > 1 , (3.19) 



i. e. in both limits, t < 1 and t»1, the viscosity coef- 
ficient is much less than the critical value given by Eqs. 
(3.15) or (3.16). The function (3.18) reaches its maximum 
at t — £>j -1 , and this maximum is given by 



cr / max 



Bi 



(3.20) 



Even in this, the least favorable case [when r ~ 0.5 as 
one can see from Eq. (2.4)] v is less than v CT by a fac- 
tor of 2 or so, which is enough for small-scale turbulence 
to appear. Therefore the range of physical conditions un- 
der which the viscosity due to fully developed turbulence 
should dominate is indeed very broad. 

Unlike the classical example of gravitational instabil- 
ity, our mechanism for the growth of shortwave pertur- 
bations has a much higher level of viscous stabilization, 
as it follows from the value of critical viscosity calculated 
above in comparison with that for a self-gravitating disk. 
Indeed, the growth of shortwave perturbations leads to the 
increase of the amplitude by a factor of 

Vi r {U) .kr(0) 



MO) 



This strong inequality follows from the fact that 
\k r {Qi)/k v \ <~ i*/T, where i* is the characteristic time of 
the growth of perturbations and T is the period of the 
disk revolution. According to the perturbation theory im- 
plemented to examine instability, the condition t*/T >> 1 
holds (otherwise the zero-approximation of perturbation 
theory is not fulfilled: the equilibrium condition is broken 
for the time less then that of one revolution of the disk). 
As a result of the above inequality, a strong growth of 
perturbation takes place, which leads to the development 
of short-scale turbulence and the appearance of turbulent 
viscosity. A very large factor of the growth given above 
explains why the value of the critical turbulent viscos- 
ity, which stops the growth of perturbations, turns out 
to be much larger than that for the instability of a self- 
gravitating disk (e.g., Fridman & Polyachenko 1984, p. 
41). 

4. Viscosity in the Circumnuclear Ring 

Before making numerical estimates, we list the basic pa- 
rameters of the clumpy gas in the circumnuclear ring 
(CNR), such as the inferred clump size a, the volume fill- 
ing factor F, and velocity dispersion of the clumps a v , 
taken from Jackson et al. (1993) and Glisten at al. (1987): 

a~0.15pc; F ~ 0.1 - 0.3; a v ~ 20 km s _1 . (4.1) 

Adopting the average gas density in the clumps n = 
10 5 cm~ 3 one finds the average clump mass m — 5 M©. 
This gives the ratio ac/a ~ 10~ 3 , which implies that elas- 
tic (gravitational) interactions between the clumps are 
negligible compared to inelastic ones, i. e. clump-clump 
collisions. In other words, gravitation plays no role in the 
interactions between the CNR clumps. 

The mean free path of the clumps given by 



4 a .„ „ 

~ 0.7 

3F v 



2) pc , 



(4.2) 



is rather large (even in a marginal conflict with the sim- 
plifying assumption (3.1) that I -C h). The clump-clump 
collision rate in the CNR is given by: 



1) • 10 



12 s- 1 



(4.3) 



i.e. uj c <fl ~ 2 • 10~ 12 s _1 implying less than one col- 
lision per revolution. The anticipated optical depth is 
t ~ 0.1 — 0.5, which, according to Eq. (3.17), results 
in v < v CT . Therefore, the conditions for fully developed 
turbulence to appear, which are described at the end of 
Sec. 3, are met to yield the viscosity coefficient 



Vturh — 1 



10 24 cmV 



(4.4) 



It is instructive to compare this result with the upper 
limit to viscosity derived by von Linden et al. (1993a, b). 
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Their results were obtained by fitting an accretion disk 
model with arbitrary viscosity to the velocity fields of vari- 
ous molecular clouds, and then inferring the required kine- 
matic viscosity from the fit. Successful fits were made in 
the radial range from 10 to 100 pc, with an implied kine- 
matic viscosity of 6 1 26 cm 2 sec -1 at a distance of 100 pc. 
Such a high viscosity is just within the limits imposed by 
the basic assumption of an accretion disk: The thin disk 
assumption implies, in the context of isotropic turbulence, 
that the kinematic viscosity has to be clearly less than the 
circular velocity times the scale height. At r ~ 100 pc the 
circular velocity is ~ 200 km/sec, and the scale height is 
difficult to determine; the z-distribution of clouds gives a 
lower limit to the scale height, and that is <~ 10 pc. This 
means that the kinematic viscosity has to be less than 
6 1 26 cm 2 sec -1 ; the fit by Linden et al. is obviously just 
at the limit. It follows either, a) that the real scale height 
is quite a bit larger, with a rather hard limit at roughly 
1/3 of the radius, implying a hard limit of the viscosity 
of 2 10 27 cm 2 sec -1 , which should not be reached, or b) 
that the kinematic viscosity cannot be described with an 
isotropic turbulence. 

The viscosity values implied by the fit to the obser- 
vations of molecular clouds differ for different radii. A fit 
was made for clouds at 100 pc as well as 10 pc, with the 
viscosity decreasing for smaller radii. This radial varia- 
tion may become steeper at smaller radii. The hard upper 
limit mentioned above would imply that it decreases as ap- 
proximately r 113 , and so would imply that the hard upper 
limit at 1.5 pc is w 5 1 24 cm 2 sec -1 . We note that the re- 
sult (4.4) is a factor of 5 below this upper limit. Assuming 
the same scaling for the kinematic viscosity derived by von 
Linden et al. would lead to w 1.5 1 24 cm 2 sec -1 , which is 
close to Eq. (4.4); assuming the number derived from a fit 
at 10 pc we obtain an estimate which is near the limit. 

5. Discussion 



A qualitative argument presented at the end of Sec. 2 
shows that the viscosity in a cloudy disk is, in general, 
much larger than that for a continuous disk. This assertion 
could be easily confirmed by numerical estimates when 
we address the circum- nuclear ring (CNR): in a continu- 
ous disk with similar global parameters, one would have 
v ~ c s l<c s h ~ 10 23 cm 2 s -1 , where c s ~ 1 km/s is the 
sound speed of the gas (Glisten et al. 1987) and h ~ 0.5 pc 
is the thickness of the disk. What provides a much larger 
viscosity than it would be possible in a continuous disk is 
(although only partly) the clump-clump collisions. As is 
shown in Sec. 3, the collective mode of instability would 
dominate the dynamics of a clumpy disk if viscosity in the 
latter is less than some critical viscosity [Eq. (3.11)] whose 
value is given by Eqs. (3.15) or (3.16). This condition is 
(marginally) met in the CNR, which results in the turbu- 
lent viscosity coefficient f tU rb — 1 10 24 cm 2 s -1 . We note in 



passing that Giisten et al. (1987) suspected the existence 
of turbulence in the clumpy CNR, although they did not 
evaluate its viscosity coefficient. 

It is worth mentioning that, in spite of the differences 
in the specific physical parameters of the clumps in the 
CNR and active galactic nuclei, or AGN (for the latter, 
see Netzer 1990), the value of t'turb evaluated above turns 
out to be in the range of the values estimated by Fridman 
& Ozernoy (1992) for cloudy disks in AGN. This might be 
relevant to the issue whether the CNR could be considered 
as a prototype for circum-nuclear tori around some types 
of AGN. (It is interesting that a large value of the viscosity 
in the CNR implies its rather large scale height, which 
is a required geometry for the AGN tori). In any case, 
the angular momentum transport in the CNR seems to 
be a template while considering similar issues both for 
quiescent and active galactic nuclei. 

The analysis performed above has not accounted for 
the magnetic field in the CNR. Meanwhile several observa- 
tional techniques revealed the field strength in the clouds 
to be ~ 1 mG (for a recent review, see Genzel et al. 1994). 
A magnetic field as strong as this cannot be ignored in the 
transport of angular momentum. Conservation of angu- 
lar momentum transported by both viscous and magnetic 
stresses can be written for the CNR in which v v w const 
in the form (Ozernoy & Genzel 1998): 



M 



27T7-E 



(t,v A + ^cfffi) 



I [ ?± 

r 



1/2 



-i -1 



(5.1) 



where M is the mass inflow rate, E is the surface density 
of the CNR, £ <~ 1 is the (—B r /B v ) averaged over the 
z-coordinate, va is the Alfven velocity, v c g is an effective 
viscosity, and Ri ~ 1.5 pc is the inner radius of the ring. 
If we consider the CNR as a magnetized disk for which 
va — hfl ~ 30 km/s and v e s is given by Eq. (4.4) the 
two terms in parentheses in the r.h.s. of Eq. (5.1) turn 
out to be comparable. This implies that while evaluating 
the angular momentum transport in the CNR, both mag- 
netic and turbulent viscosity need to be accounted for. 
One can see that with the parameters listed above and 
S ~ 2 10 -2 g cm -2 the total mass inflow rate given by 
Eq. (5.1) at a fiducial distance of r = 2 pc amounts to 
M ~ 10 -2 Mq yr -1 . 

This result is consistent with a naive, by order-of- 
magnitude, estimate of M ~ M/t tr , where M is the 
CNR mass and t tr = ^f/tWb is the characteristic time 
for the angular momentum transport: Taking f tU rb — 
1 10 24 cm 2 s -1 and M <~ 3 10 4 M Q , i.e. somewhere in 
between 10 4 and 10 5 M Q , the current estimates for the 
CNR mass (Genzel et al. 1995), one gets tt r — 6 10 5 yr 
and M ~ 5 10 -2 M Q yr -1 . This is consistent as well with 
the inflow rate toward the Galactic center inferred from 
the observational data (e.g. Blitz et al. 1993, Genzel et al. 
1994). 
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Appendix A: 

Shear Viscosity in a Clumpy Disk at Elas- 
tic/Inelastic Encounters 

The shear viscosity in a differentially rotating clumpy 
disk, with accounting for gravitational interactions be- 
tween the clumps of mass m, is given by Eq. (60) in Stew- 
art & Kaula (1980): 



-ril'v = A a 2 v , 



A 



9^2 nuir 1 



(36/5) + 55 (rJw G 1 ) 2 ' 



(Al) 



where fl' = dO/dr, and loq is frequency of clump-clump 
collisions given by 



LO G = 



3.35 G 2 m 2 n 



(A2) 



Eq. (A2) corresponds to Eq. (57) of SK when one puts 
ln[l + (al/GmVl) 2 ] 1 / 2 = 1 because in a flat disk the long- 
range gravitational interactions could be neglected com- 
pared to the short-range ones. The value of log is related 
to the characteristic time of gravitational encounters, t Gl 
given by Eq. (2.3) simply by lOq 1 = ta/V2. By substitut- 
ing Eq. (A2) into Eq. (Al), the latter could be written in 
the form: 

v=4 AG \ , (A3) 
n (1.95r) 2 + l V ' 

where Aq = 0.83 for a Keplerian disk, Aq — 1.25 for a 
solid-body rotating disk, and r = flto- 

The viscosity coefficient for clump-clump collisions, 
when the clump gravity is negligible [Goldrcich & 
Tremaine 1978, Eq. (46)], can be written in the form anal- 
ogous to Eq. (A3): 



0.46 



(0.97 rf + 1 



(A4) 



Eq. (2.1) unifies the representation of viscosity both in 
the case when interactions are inelastic, while gravitation 
is negligible [Eq. (A4)] and in the case when gravitational 
encounters are dominating [Eq. (A3)]. 

Appendix B: 

Character of Perturbations in a Differentially Ro- 
tating vs. a Solid-body Rotating Disk 

Let us consider an incompressible, differentially rotat- 
ing, liquid cylinder. The linearized equations of motion 
take the form: 







divu 


= o, 


dv r 
~dt 


OLf 


- 2m v 


1 dP 

p dr ' 




+ 


K 2 „ 
+ V r 

2Vl r 


1 dP 

pr dip 



Substituting a perturbation in the form of the wave 
packet / ~ exp[— iuj + mip + j k r dr] results in the following 
dispersion relation (k v = m/r): 



i(2Q- 



n 2 /2Q)^f. 



As always k 2 < 4fl 2 , and therefore the condition of in- 
stability is k r k v < 0, which corresponds to the growth 
of the leading spirals. Under this condition the instability 
is monotonous without any oscillations with epicyclic fre- 
quency. If the rotation is solid-body, one obtains to = 0, 
i.e. the perturbations are stationary. 

The above solution has a drawback since the r- 
componcnt of the group velocity for perturbations under 
consideration is zero. A more detailed analysis should be 
performed to find out what is the role of r-dependence 
for Q. This difficulty can be overcome if one considers the 
system in a corotating frame of reference. In this case, an 
exact time-dependent solution can be derived (see Frid- 
man, 1989) which demonstrates a qualitatively the same 
behaviour as is obtained above. Specifically, in any dif- 
ferentially rotating disk whose angular velocity decreases 
with radius, the leading spirals grow monotonously. As for 
a solid-body rotating disk, all perturbations are station- 
ary. 

Note that if one would consider absolutely compress- 
ible case, i.e. P = in the above equations, the continuity 
equation takes the form: 

dp ,. _ dp 

— + pdivv + v r — = 0, 

ot dr 

and we immediately obtain 

- 2 2 
U) = K . 

The contrast to other work is illuminating. Julian & 
Toomre (1966) considered perturbations in a very thin 
collisionless self-gravitating stellar disk. In fact, for such 
perturbations a collisionless stellar disk is absolutely com- 
pressible, in contrast to incompressible character of per- 
turbations with wavelengths much smaller than the disk's 
width for a gaseous disk considered here. Goldreich & 
Lynden-Bell (1965) studied dynamics of perturbations in 
a gaseous disk, but only at high frequencies. Really, per- 
turbations considered in their paper demonstrate (asymp- 
toticaly) oscillations with sound frequency kc, where k is 
the wave number and c is the sound speed. As the scale of 
perturbations is smaller than the disk width, i.e. kh > 1, 
their frequency is higher than the epicyclic one. But as 
the equations of motions contain three derivatives with 
respect to time, they describe three modes of perturba- 
tions. Besides the two sonic high frequency modes, one 
mode of low frequency perturbations exists as well. What 
is examined in our paper is the role of low frequency mode 
in dynamics of the circumnuclear disk. As was shown in 
Fridman (1989), this mode does not demonstrate any os- 
cillations and is stationary in a solid-body rotating disk. 
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Appendix C: 

Velocity Perturbations in a Plane Shear Layer 

We start from linearized equations for short-wave, low- 
frequency perturbations in a barotropic gaseous disk in the 
corotating frame of reference (for derivation, see Fridman 
1989): 



dVr 00 - 



0, 



-ik r (t)x — vk\(t)v r 



+ (2il -A)v r = -ikpic ~ vk 2 ± (t)v, 



v 



Here 51 = const is the angular velocity of the corotating 
frame of reference; A = —r a (dVL/dr) a characterizes the 
value of the shear; k r (t) = k r + tAk v ; k 2 ± (t) = k 2 (t) + 

This system has an exact solution in the form (Frid- 
man, 1989): 



v r (t) = 



k\{t) 



v r (0) exp 



t 

o 



j_{t)dt 



(CI) 



For small t the exponential term is negligent and we come 
up with vorticity conservation: 

curl z (v) = i(k r v ip — k ip v r ) = — -^-k"]_v r = const. (C2) 

For a plane shear layer, a similar set of equations can 
be derived. The simpliest way to do this would be to put 
in the above system f2 — > 0, r — > oo, fir — > v(x). As a 
result, we obtain the following system: 



k x {t)v x ~t~ kyVy 

dv x 
dt 



0, 



dt 



Bv x = 



-ik x (t)x - vkj_(t)v x , 
-ikyX- vk\(t)v y . 



Here B = (dv/dx)o characterizes the value of the shear; 
k x (t) = k x — tBk y ; fc^(t) = k 2 (t) + k 2 . One can easily 
see that this system is almost similar to the previous one 
and it has a similar solution: 



v x (t) 



k 2 ± (t) 



v x (0) cxp 



-v J k 2 ± (t)dt 



(C3) 



Obviously, in a plane layer there are no any Coriolis forces 
at all. Nevertheless, the solution (C3) has the same form 
as (CI). What is relevant to the real physical meaning of 
our solution is conservation of vorticity (eq. C2), which 
has been emphasized here. 
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sented by Eq. (3.3) for several values of (i = fc r (0)/fc v . The 
solid lines show a non-viscous case (Eq. 3.6). A general, 
non-zero viscosity solution (3.3) is shown by the dotted 
lines for v = 0.3 f cr it and by the dashed lines for v — ^ cr it, 
where ^ cr ;t is defined by Eq. (3.5). 



References 

Bcgelman, M.C, Frank, J., and Shlosman, I. 1989, in "Theory 
of Accretion Disks" (eds. Meyer, F., Duschl, W.J., Frank, 
J., & Meyer-Hofmeister, E.). Kluwer, Dordrecht, p. 373 



